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Abstract 

We compute the one-loop radiative corrections to the masses of the top quark, the 
stop squarks and the gluino in the Minimal Supersymmetric Standard Model. We in- 
clude the loops controlled by the strong coupling constant and the top Yukawa coupling, 
neglecting those controlled by the SU{2) x U(l) gauge couplings and by other Yukawa 
couplings. We find that the significant scale-dependence of the renormalization- group- 
improved tree-level expressions is almost completely removed. Even for natural choices 
of the renormalization scale, corrections can be numerically relevant. Our results should 
allow more reliable predictions to be extracted from candidate fundamental theories 
beyond the MSSM. 
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1 Introduction 



The evaluation of radiative corrections in the Minimal Supersymmetric Standard Model 
(MSSM) plays an important role in our understanding of its phenomenological implications. 
This is of particular interest, the MSSM being, at the moment, the most promising extension 
of the Standard Model (SM) for the strong and electroweak interactions. 

Finite corrections to supersymmetric Ward identities, such as those associated with 
Higgs boson masses and couplings 0, have a direct impact on the model predictions. The 
remaining corrections have more subtle phenomenological implications, since they can be 
largely reabsorbed into redefinition of the many arbitrary parameters of the MSSM. This 
is no longer true, however, if one embeds the model into a more fundamental theory (as 
superGUTs, super gravities or superstrings) where some or all of the MSSM parameters can 
be predicted. 

In this paper we compute the one-loop corrections to the masses of the top quark, the 
stop squarks and the gluino, going beyond the leading-log approximation. We present ex- 
plicit analytical formulae, keeping in the Feynman rules only the vertices proportional to 
the strong coupling constant g s and to the top Yukawa coupling h t . This approximation 
should reproduce the main contributions to the one-loop radiative corrections. Comparing 
our full one-loop results with the renormalization-group-improved tree-level masses, we find 
that there can be significant differences, also in the case of "natural" choices for the renor- 
malization scale. The corrections can be particularly relevant for the lightest stop eigenstate 
and for a very light gluino. In the top quark case, we find that the extra corrections due to 
loops of supersymmetric particles can be non-negligible compared to the pure QCD ones. 

The paper is organized as follows: in this section, we introduce the notation by giving the 
top, stop and gluino tree-level mass formulae; in section ^| we introduce the renormalization 
procedure, we recall the relation between running and pole masses, and we give the main 
points of our calculations; in section |3] we discuss the numerical results and their phenomeno- 
logical implications. The complete results for the one-loop-corrected top, stop and gluino 
masses are presented in appendix A. The explicit solutions of the one-loop renormalization 
group equations for all the parameters of relevance to our calculations, in the approximation 
where the only non-zero dimensionless couplings are g s and ht, are reported in appendix B. 

The MSSM [0 contains three generations of quark (Q { , Uf, Df) and lepton (L i; Ef) chiral 
superfields (i = 1, 2, 3), and two Higgs chiral superfields (Hi, H 2 ), with interactions specified 
by an i?-parity-invariant superpotential and a collection of soft supersymmetry-breaking 
terms. As often done in the literature, we shall assume real mass parameters and neglect all 
Yukawa couplings apart from the top quark one. The tree-level expressions for the different 
masses and couplings in the MSSM are well known: in the following, we shall use the notation 
and conventions of ||, including the reversal of the sign of [i in the squark sector as specified 
in the Errata ||. 

We now recall the tree-level expressions for the top and stop masses (the tree-level gluino 
mass is just the explicit soft mass itself). With no loss of generality, the Vacuum Expectation 
Values (VEVs) of the Higgs fields, i>j = (H®) (i = 1,2), can be chosen to be real and positive. 
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The top quark mass, then, is given by: 



1,2 2 

h t v 2 . 



The mass matrix for the stop squarks is: 



m t + m Q + (l m w ~ l m z) cos 2(3 
m t (A t — p, cot (3) 



in 



m t (A t - ii cot (3) 
2 1 fh\ + §(m,| — m 2 ^) cos 2/3 



(1) 



(2) 



where rfiQ, rfiT, A t are soft mass terms, tan /3 — v 2 /vx, \i is the superpotential mass parameter 
and, at the tree level, 



rn 



The eigenvalues of @ are: 

1 



9 / 2 , 2\ 



1 



m z = («! +u 2 ). 



(3) 



± 



^( m Q - mr) + ^(8m^ - 5m|) cos 2(3 



+ rrit(A t — fj, cot (3)' 



The corresponding eigenstates are: 

ti = cos + sin <f) t t R , t 2 = - sin (b t t L + cos t t fl , 
where the rotation angle <p t that diagonalizes the mass matrix is given by: 

m t (A t — fi cot (3) 



tan! 



[\{m 2 Q - fh%) + ^{Sm^ - 5m|) cos 



(4) 



(5) 



(6) 



with — ~ < <f> t < |. Our conventional ordering of the eigenvalues (namely, m? > m 2 2 ) 
is obtained with the prescription that 4> t should have the same sign as the combination 
A t — /icot (3. 

Summarizing, at the tree level the top, stop and gluino masses are functions of the 
following independent input parameters: hf, g, g , ,mz,m g ,mQ,fhr, A t , // and tan/3. As we 
shall see in the next section, other parameters will appear via the one-loop corrections. 



2 One-loop renormalization 

At the one-loop level, the pole masses of the top and the gluino are given by: 

m t = [m t {Q)-UMQ)y =mt , (7) 

™ a = [m 9 (Q)-^-MQ)]\j/=fn a - ( 8 ) 



2 




Figure 1: One-loop self-energy diagrams for the top: a) iHf; b) iUj; c) iHf ; d) illf ; e) 

mf. 

In the above equations, Q is the renormalization scale; mt,g(Q) are the running masses in 
some renormalization scheme, while Ut, g (j/,Q) are the real parts of the renormalized self- 
energies in the same scheme. Up to higher-order corrections, the pole masses do not depend 
on the scale or on the renormalization scheme. 

To evaluate the corrections to the pole masses of the stops, we need to study the corrected 
inverse propagator matrix T^. In the (ti,tz) basis, where ti are the tree-level eigenstates, 
the expression for the inverse propagator is: 



y ( p 2 -<(g) + n n (p 2 ,Q) n 12 ( P 2 ,Q) 
13 \ n 21 ( P 2 ,g) p 2 -m 2 2 (Q) + n 22 ( P 2 ,Q) 

At the one loop, the stop pole masses are: 



(9) 



<=<(Q)-nn(p\Q) „,_^ 2 , (10) 

™2 ~ 2 



m 2 t =mi(Q)-U 22 (p 2 ,Q) 



• 2 ' 

-mi 



*2 
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In eqs. (^), (j|), ( p!0|) and (PD, the running masses can be calculated using the tree-level 
functional forms (presented in section [I]), with all the parameters considered as Q-dependent. 

To evaluate the one-loop self-energies of top, stops and gluino, we work in the approxi- 
mation where all dimensionless couplings other than g s and h t are set to zero, which should 
take into account the most relevant one-loop corrections to the studied observables. The 
Feynman graphs that represent the self-energies of top, stops and gluino in this approxima- 
tion are reported in figs. [I], @ an d ^Q- The self-energies for the top (Hit) and gluino (iHg) 
read: 

iIL t (jfi = iUl + iUl + iUf + mf +iUf + , (12) 

m s (tf) = /ii;; • ;u' :l ■ /ii!; . (13) 



1 The graphs in fig. |^b and ||c must be counted twice, due to the fact that the gluino is a Majorana 
fermion. 



3 




Higgses 



stop 



Figure 2: One-loop self-energy diagrams for the stops = 1, 2): a) iHfj; b) ill*,-; c) zllf • 



while for the stops the self-energies iHij can be written as: 

fflnGo 2 ) = ini 1 + m[ 1 + m\ 1 + Mf 1 + m& + iiig, (u) 

m 22 ( P 2 ) = in| 2 + «n' 22 + m b 22 + mg + mf 2 + TTf 2 + , (15) 

ffl 12 ( P 2 ) = ziif 2 + inf 2 + m5 2 + *nf 2 + ing + m* 2 , (16) 

ffl 21 ( P 2 ) = tp 12 ( p 2 ). (17) 

They have been grouped in subsets of graphs, to make the results more transparent and 
to simplify the numerical study that will be presented in the next section. The masses 
of the particles circulating in each subset are controlled mainly by one input parameter, 
e.g. the contributions of the graphs with Higgs and Higgs/squark loops to the stop self- 
energies (fig. ^]d) can be studied by varying m\, while neutralino and chargino loops (figs, 
^e and 0f ) depend mainly on \i. The self-energies have been calculated both in the Feynman 
and in the Landau gauge, with the same results. As usual, in order to explicitly preserve 
supersymmetry, we have adopted the DR scheme ||. The complete results for the one- loop 
self-energies of the top, the two stops and the gluino are reported in appendix A. 

It can be seen by looking at the formulae that some new parameters must be introduced 
in order to calculate the full one-loop expressions. The independent parameters that must 
be added to the ones that appear in the tree-level relations are a s , and %, where is 
the physical mass of the CP-odd neutral Higgs boson and ffiB is the soft mass term for the 
right sbottom squark. The masses of the other Higgs bosons (H°, h° and H^) and the CP- 
even mixing angle a can be deduced with tree-level relations from mA,mz,mw and tan/3. 
The masses of the sbottoms can be obtained from mQ,mB,mz,m\v and tan/3. All these 
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gluino 



Figure 3: One-loop self-energy diagrams for the gluino: a) ill-; b) zTl|; c) ill- 



tree-level relations can be found in the literature. 

We have verified that the inclusion of terms proportional to g and g' (namely, m 2 z and 
rriyv) in the tree-level relations for the masses of the particles circulating in the loops does 
not significantly modify the results. For this reason, we have neglected them in the tree- 
level relations for the charginos and neutralinos mass matrices. This further approximation 
greatly simplifies the numerical implementation of the general formulae reported in appendix 
A. The charginos and neutralinos eigenvalues in this approximation are: 



m xt 



(M X ,M 2 , n,-y) j 



1,2, 
1,4, 



(18) 
(19) 



while the rotation matrices are: 



Vi 



Ni 







1 

~V2 



(20) 



(21) 



V2 I 



In this approximation it can be seen, by looking at the formulae reported in appendix A, 
that the gaugino soft masses Mi, M 2 do not appear explicitly. Therefore they have not been 
included in the input parameters set. 

In our approximation the electroweak gauge couplings g and g' (equivalent to sin 2 9w and 
a em ) do not evolve under the renormalization group equations. They have been fixed at the 
values sin 2 9w = 0.23 and a em = Another fixed input parameter is the Z° physical mass, 
irtz = 91.18 GeV. The other two couplings (a s , h t ) have been chosen in order to give, at the 
electroweak scale, a value for the top pole mass that lies approximately in the experimental 
range. We stress that all the input parameters must be interpreted as MSSM parameters: 
e.g. a s is the MSSM strong coupling constant, and not the SM one. 

In order to compare our full one-loop calculations with the results obtained in the leading- 
log approximation via the renormalization group equations, we have used the explicit solution 
of the RGEs that can be found in their general form in ||. In our approximation we must 
neglect terms proportional to dimensionless couplings other than g s , h t also in the RGEs, 
which therefore can be explicitly solved. The details of the RGEs solutions are reported in 
appendix B. It must be noted that in those formulae another two parameters (deducible from 
the already discussed set) appear, v 2 and m# 2 . For consistency, their expression in terms of 
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Figure 4: One-loop self-energy diagrams for the Z° : a) —ig^H^; b) —ig^U^; c) —ig^lV^ ; 
d) -ig^Tlf . 



the input parameters has been calculated at the one-loop level. To extract from the input 
parameters the one-loop VEV V2, we have calculated the running DR mass m 2 z (Q) from the 
physical Z° mass mz using the following relation: 



m%{Q) =m|+ Tl z (p 2 ,Q) 
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[22) 



In the already obtained expression for the one- loop self-energy Hz 0, we have kept only 
terms 0(g 2 ) to implement our approximation^. The Feynman graphs that contribute to the 
Z° self-energy in this approximation are reported in fig. |j. The decomposition of Hz on the 
subsets of graphs is: 

- ig^iiziv 2 ) = -ig^ [n z ( P 2 ) + nfV) + n§V) + nf ( P 2 )] . (23) 

The one-loop expression is presented for completeness in appendix A, jointly with our new 
one-loop results. The one-loop VEV t>2 is then obtained from m 2 z (Q) and tan {3 with the tree- 
level relation (Q). This procedure corresponds to the minimization of the one- loop effective 
potential: with this definition of the VEVs, the one-loop tadpole terms disappear (as the 
tree-level tadpoles disappear as a consequence of the tree-level scalar potential minimization). 
To extract from the input parameters the one- loop value for m^ 2 , we have used the one- loop 
expressions that relate mH 1 ,mH 2 and mz,rriA presented in eqs. (22), (23) and (24) of |§. 

In the case of the stop squarks we have verified, as a check of the calculation, that 
the quadratic divergence of the sum of the full set of graphs is zero, as it should be in a 
softly broken supersymmetric theory. As a further check, we have verified analytically for 
all the studied masses that the implicit Q-dependence of the running mass cancels with the 
explicit Q-dependence of the one-loop self-energy. The logarithmic dependence on Q for the 
renormalized gluino mass has also been checked by a comparison with the results in . The 
one-loop correction to the top quark mass (including the finite part) in the DR scheme for 



the terms proportional to g s has been compared with |]IU||. As a consistency check, we have 



2 Of course it is not possible to neglect all the vertices proportional to g and </, for this will result in 
putting U z = 0. 
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verified that if the off-diagonal contributions 0(h 2 ) in eq. ([]) are kept the one-loop-corrected 
pole masses do not significantly change. 

3 Numerical analysis 

In this section we describe in some detail the numerical results obtained for the one-loop- 
corrected masses of the top, the two stops and the gluino. 

In fig. [| we report the scale-dependence of the DR running masses (dashed lines) and 
of the one- loop-corrected pole masses (solid lines). Our representative choice of parameters, 
assigned conventionally at the scale Qo = 1000 GeV, is: rriQ = rhx — rhs — rh g = A t = fi = 
mA = 200 GeV; tan/5 = 2; it can be seen that, in all four cases, the strong Q-dependence 
of the running mass m(Q) is significantly reduced. This is an expected feature, due to the 
analytic cancellation of the logarithmic dependence in Q between the running expression 
and the one-loop self-energy. It must be noted that the residual scale-dependence present in 
figs. ||b and |^c can be explained in the following way: using in the numerical calculations the 
neutralino and chargino masses reported in eqs. (|18|) and (|T9|), we introduce a partial non- 
cancellation of the scale-dependence between fh^.(Q) and Ylu^rhf., Q)f\ Due to the improved 
scale-independence, one can check if the usual procedure adopted to choose the optimal scale 
to evaluate the running mass (Q ~ m) can lead to a wrong calculation of the one-loop effects. 
This is actually the case, as can be seen for example for the gluino or the lightest stop mass. 

Now we can start a detailed analysis of the complete results, comparing them with the 
running masses obtained at the "natural" scale Q ~ m. We find that in many cases the 
largest corrections to the running masses come from the loops in which a gluino is circulating. 
Therefore we report, for all four masses under consideration, the dependence of the one- loop- 
corrected mass on the gluino mass. We present also the dependence on other parameters in 
various interesting cases. 

1. In fig. ^| we report the one-loop-corrected mass for the top. We have verified that 
the top mass corrections are almost independent from m^, which has then been kept 
fixed at m>A = 300 GeV. The corrections have been calculated at the "natural" scale 
of Q = 170 GeV. We have studied two main situations: the case in which the two 
stop states are almost degenerate and the case in which the lightest stop state is 
"very" light (between the experimental lower bounds of ~ 50 GeV and 100 GeV). To 
obtain degenerate stop states, we have chosen universal soft mass terms (rhq = my = 
rfiB = rhu) for the squark and A t = fi cot (3 = 250 GeV. This case is illustrated in 
fig. |a, for four representative values of rhu, ('rhu = 100,200,400,600 GeV). In the 
same situation, we have studied the dependence on tan (3. The results are reported 
in fig. |6|b, with fh g = 300 GeV, tan/3 =1-50. To obtain a light stop state, we have 
still chosen universal soft mass terms (mg = my = mg = rhu) f° r the squark, but 

3 By studying the Q-dependence of the different sets of Feynman graphs that appear in fig. @, we have 
seen that the scale-dependence is actually reintroduced adding graphs (e) and (f) to the others. 
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Figure 5: Scale- dependence for the one-loop-corrected masses: a) t; b) t\; c) t 2 ; d) g. The 
solid lines represent the complete one-loop results, while the dashed lines are the running DR 
masses. 



A t = 0. The dependence of the top mass on itla in this situation is reported in fig. 
|c, with m A =200-600 GeV. In these last two figs. m v = 200, 400, 600 GeV. It can 
be seen that in all the cases presented, the corrections to the running top mass are 
quite stable and not very sensitive to the variation of the parameters. The results 
obtained for rhjj = 150 GeV are reported in fig. |5|d. The corresponding lightest stop 
mass is approximately m t , 2 < 100 GeV. In this last figure, we have studied in detail the 
single contribution of the different classes of Feynman graphs to the complete result, 
in the case of a "very light" t 2 . The corrections have been compared with the QCD 
correction (represented by the solid line that lies in the middle), which is the main SM 
contribution in the approximation g — g' — 0. The QCD correction to the running 
mass at Q — 170 GeV is slightly less than 10 GeV. The dotted lines represent the 
different contributions added to the QCD corrections. Starting from the QCD line, we 
can see the charginos, neutralinos, Higgses and gluino contributions. It can be seen 
that the largest correction (to the QCD result) comes from the gluino loopsQ. We stress 
that, in the different cases examined, the MSSM corrections can reach half (or more) of 
the QCD correction. Also, the complete one-loop results can significantly differ from 
the running ones. 

2. In fig. [7| we report the results for the "heaviest" stop state i\. Also for this particle 
we have studied the two cases of almost degenerate stop states and of a "very light" 
t 2 . In fig. 0a we consider the case of degenerate stop states at three different values 

4 Stop/gluino loops have already been taken into account in Jicf] , with similar results. 
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Figure 6: One-loop-corrected top mass (in GeV): a) as a function of the running gluino mass, 
for m Q = m T = m B = 100, 200, 400, 600 GeV from the bottom; A t = 250, /i = 500, m A = 300 
GeV; tan/3 = 2. b) as a function o/tan/? ; for rhq = rh T = rh B = 200,400,600 GeV 
from the bottom; A t = 250, /i = 500, = 300; fh g = 300 GeV. c) as a function of tua, 
for rfiQ = m T = m B = 200,400,600 GeV from the bottom; A t = 0,/j, = 500, m g = 300 
GeV; tan/3 = 2. d) as a function of the running gluino mass, for rriQ = friT = rh B = 
150 GeV; A t = 0, fi = 500, = 300 GeV; tan/3 = 2. The solid lines are the complete 
result, the dashed ones represent the running mass. In (d), the lower solid line is the pure 
QCD correction; the dotted lines represent the contributions of (starting from the bottom) 
charginos, neutralinos, Higgses and gluino loops added to the pure QCD correction. 



of the universal soft mass term fhjj = 100, 200, 400 GeV. As usual the dashed lines 
represent the running masses. The three results have been obtained at the "natural" 
scales Q = 150,240,430 GeV. In figs. 03 and |7|d we have reported on the dependence 
on rriA and /i in the case of a "very light" t 2 at the scale Q = 300 GeV (with rhjj = 170 
GeV, A t = 0,m g = 300 GeV). In particular, in 0b fj, = 500 GeV and m A =200-600 
GeV while in |7|d = 500 GeV and fi =200-600 GeV. In fig. 0c, is reported on the 
dependence of rh tl on the gluino mass in the same situation, with fhu = 200 GeV 
and Q = 320 GeV. In this case, we have studied also the contributions coming from 
the different loops. The horizontal dotted lines are (from the bottom) the neutralinos, 
charginos, Higgses, sbottoms and stops contributions, while the dotted line that follows 
the complete result is the gluino contribution. It can be easily seen that also in this 
case the most important contribution comes from the gluino loops. Actually, in all 
the cases that we have observed, the corrections to the m tl running mass are not very 
large. We have seen that for almost degenerate running states, the one-loop-corrected 
masses can be interchanged (namely, i 2 can become the heaviest stop state). 
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Figure 7: One-loop-corrected ti mass (in GeV): a) as a function of the running gluino mass, 
for m Q = m T = m B = 100, 200, 400 Ge V from the bottom; A t = 250, fi = 500, m A = 300 
GeV; tan/3 = 2. b) as a function of niA, for rhq = rh T = rh B = 170 GeV; A t = 0,/i = 500, 
frig = 300 GeV; tan/5 = 2. c) as a function of the running gluino mass, for rfiQ = rfiT = 
friB = 200 GeV; A t = 0, /i = 500, = 500 GeV; tan/3 = 2. d) as a function of \i, 
for ffiQ = friT = friB = 170 GeV; A t = 0,niA = 500, m g = 300 GeV; tan/3 = 2. The 
solid lines are the complete result, the dashed ones represent the running mass. In fig. (c) 
the horizontal dotted lines represent the contributions of (starting from the bottom) Higgses, 
sbottoms, charginos, neutralinos and stops loops. The dotted line with the same shape as the 
complete result is the gluino contribution. 



3. In fig. H we report the results for m t2 with the same choices of parameters as for 
fh tl . In fig. |8|a, it can be seen that £2 has a behaviour quite similar to t\\ actually, 
the one-loop corrections can interchange the two states (£2 can become the heaviest 
one). Also for £2, the corrections to the running mass (evaluated at the same "natural" 
scales) are quite small. In figs. ||b, ||c and |]d we have studied the case of a "very 
light" t 2 , at the "natural" scale Q = 100 GeV. It can be seen in all three figures that 
the one-loop corrections significantly modify the renormalization-group-improved tree- 
level results. The most interesting situation is reported in fig. [5]c. At the "natural" 
scale of Q = 100 GeV, we have chosen a universal soft mass term fhy = 200 GeV, in 
order to get a sufficiently heavy running mass, fht 2 (Q) — 150 GeV. It can be seen from 
the figure that the one-loop-corrected mass is significantly different from the running 
mass, and that can easily fall below the experimental lower bound of ~ 50 GeV also 
for a non-super-heavy gluino mass. This implies that in the case of a light stop state 
the running-mass result must be used (also at the "natural" scale) with great care. 
As usual, we have studied the different contributions to this case. Starting from the 
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Figure 8: One-loop-corrected i% mass (in GeV): a) as a function of the running gluino mass, 
for m Q = m T = m B = 100, 200, 400 Ge V from the bottom; A t = 250, fi = 500, m A = 300 
GeV; tan/3 = 2. b) as a function of niA, for rhq = rh T = rh B = 170 GeV; A t = 0,/i = 500, 
frig = 300 GeV; tan/5 = 2. c) as a function of the running gluino mass, for rfiQ = rfiT = 
friB = 200 GeV; A t = 0, fi = 500, = 500 GeV; tan/3 = 2. d) as a function of fi, for 
friQ = friT = ffiB = 170 GeV; A t = 0,rriA = 500, m g = 300 GeV; tan/3 = 2. The solid 
lines are the complete result, the dashed ones represent the running mass. In fig. (c), the 
horizontal dotted lines represent the contributions of (starting from the bottom) neutralinos, 
charginos, sbottoms, stops and Higgses loops. The dotted line with the same shape as the 
complete result is the gluino contribution. 



bottom, the horizontal dotted lines are the contributions coming from neutralinos, 
charginos, sbottoms, stops and Higgses loops. The dotted line with the same shape as 
the complete result represents the gluino contribution. 

4. In fig. ^| we report the results for the one-loop-corrected gluino mass in two different 
casesi a heavy gluino and a massless gluino. Our results agree with the one-loop 
corrections calculated in [O, and contain as a particular case the results of [[HJ for a 



gluino massless at tree level. In |9]a and |9|b we have studied the dependence on /i and 
tan j3 of the corrections to a heavy gluino mass, for a running mass fh g = 300 GeV at 
the "natural" scale Q = 300 GeV. It has been found that the corrections increase with 
increasing scalar soft mass terms rhQ,fhT,friB, and with decreasing tan/3. They can 
be of order 10-30% of the running mass, for stop masses not exceeding 1 TeV. In figs. 
||c and ||d we report results for the massless case (for which we have used a s = 0.15). 
Here the situation is reversed, and the corrections increase with decreasing scalar soft 
mass terms. It is well known that for a massless gluino the corrections at one loop are 
proportional to the mass differences between the squarks. This has been verified by 
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Figure 9: One-loop- corrected gluino mass (in GeV): a) as a function of fi, for nig = = 
m B = 200,300,400 GeV from the bottom; A t = 0,m A = 500 GeV; tan/3 = 2. b) as a 
function of tan f3, for rhq = rh T = m B = 200, 300,400 GeV from the bottom; A t = 0,niA = 
500, ft = 500 GeV. c) as a function of fi, for rriQ = rhx = rhs = 300,200,150,100,50 
GeV from the bottom; A t = 0,mA = 500 GeV; tan/3 = 2. d) as a function o/tan/3, for 
m Q = m T = m B = 300,200, 150, 100 GeV from the bottom; A t = 0, m A = 500, fi = 500 GeV. 
The solid lines are the complete result, the dashed ones represent the running mass. 



giving non-universal soft mass terms, keeping the average scale my at a small value. 
We have seen that, at the "natural" scale Q = 1 GeV, the corrections to a zero mass 
gluino can reach 3 or 4 GeV. Comparing our results with those presented in |12|], we 
note that we have explored a different region of the parameters spaceQ. 
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5 As long as we take non-zero gaugino tree- level masses Mj and M 2 , we can avoid the experimental bounds 
induced on /i by the non-observation at LEP of charginos lighter than 50 GeV that were considered in Oj . 
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Appendix A: One-loop results 



In this appendix we report all the explicit results for the one-loop radiative corrections to the 
pole masses of the stops, of the gluino and of the top quark, obtained in the approximation 
where all the dimensionless couplings other than g s and ht are set to zero. The corresponding 
Feynman rules have been deduced from the rules reported in || || . It must be noticed that 
(in the last of 0) a factor | is missing in the vertices between two neutral Higgses and 
two squarks. This factor, due to a combinatorial factor of 2, cancels if the external legs 
are the two neutral Higgses giving the Feynman rule reported in the same reference. The 
same happens to the two-gluon/two-squark vertex. In the formulae we have made use of the 
following integrals: 



F(l,2,3 
#(1,2,3 
G(l,2,3 
#(1,2,3 
L(l,2,3 
M(l,2,3 
iV(l,2,3 



Q 2 



Tie £ dz[m\(l -z)+ m\z\ log f^r^) 
Ke £ dz[mlz{l - z)]\og (^^) 
TZe J dz[m\z(l — z) — m\(l — z) — m\z\ log ^- ^ ' ' ^ 

nef Q dz[ m l(z-2f]\o g (p^^ 

f 1 (Til 2 3)\ 

Tie dz[m 3 z]log ( — ^ — 1 

He £ dz[m 3 (l - z)]log ( T{1 Q l ,3) ) 



where T(l, 2, 3) = m\(l — z) + m\z — m\z(l — z) — it 1 . The analytical expressions for these 
integrals can be found following [f[3|. They all depend on: 



F(l, 2, 3) = -2 + log f 0%*) + <<^,o 3 (151) 

\ Q z ) mi \m 2 J 

H — -o\f\ (mi + m 2 ) 2 - ml II (m 1 + m 2 ) 2 - ml 1/(1, 2, 3) 



where 



/(I, 2, 3) 



log 



\J (m\+m2) 2 —mj i ~ {m\—m-£) 2 —m 



yj (mi+rri2) 2 — m 2 + y/ (m\—m2) 2 —m 



2arctan ^-(^-^)^ 



log 



yj m^ — (mi—m,2) 2 + yj m^ — {mi+mz) 2 
\J in 2 — {rri-i~rri2) 2 — \J m 2 — (rnx^m^) 2 



if m\ < (mi — 777,2) 



if (mi — m 2 ) 2 < ml < (mi + m 2 ) 2 (24) 



if 777n > (mi + m 2) 
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The other integrals have the following analytical expressions: 



#(1,2,3) 



2 2 

mf — m 2 

ml 



2 / , m l 



m log — - — 1 



Q 2 



+ 



2 i 2 

7?^ + m 2 



G(l,2,3) 

#(1,2,3) 
£(1,2,3) 

M(l,2,3) 
iV(l,2,3) 



m 



1 



18 6 



'm| — m 2 \ 2 



ml log ^ + m\ log 



F(l,2,3) 



mj 

Q 2 



+ 3^(1,2,3) 



- -[m? + m^- 7^(1, 2, 3) 

= G(l,2,3)-if(l,2,3) 
/ 3m? x 



4 



2 2 

mf — m 2 



3ml 



I F(l, 2, 3) - , 
mf — m 2 / mf — m 2 



tf(l,2,3)-G(l,2,3) 



F(l,2,3) 



m 3 



2 2 

mf — m 2 



#(1,2,3) 



m 3 F(l,2,3)-M(l,2,3) 



All these integrals are symmetric under the permutation 1 
for M(l,2,3) 



2, made exception for L and 



N(2, 1,3). It has been stressed ||T3| that if some entries of the integral F are 
zero, the right expressions are not deducible by just taking the limit from these formulae. 
This happens because, if T(l, 2, 3) has some zero entries, some of the poles in the logarithm 
may disappear. Here we report the main particular cases, for completeness: 



F(l,2,0) 
F(1,1,0) 



-1 + 



1 



2 

mf 



m 2 



2 i Tfl^ o i 

m 1 log — - — m 2 log 



Q 2 2 *Q 2 



log 



Q 2 



The particular cases for the other integrals are deducible by taking the limit for the reported 
expressions and using the particular expressions for the F integral. This is not true for the 
L, M, N integrals in the case m\ = m 2 : but also in the pathological case their decomposition 
over the main integrals F, H, G is easy to obtain. Below are reported the expressions for 
the one-loop radiative corrections to the pole masses of the stops, the gluino and the top 
quark. In these expressions, N is the number of colours (N = 3), (a, (3) = 1,2,3 and 
(a, b) = 1,...,8. Moreover, C 2 (R) = N 2 ^ 1 and C 2 (G) = N. We have used a subscript c 
for the up-type quarks and squarks different from top and stop, and a subscript b for all 
the down-type quarks and squarks. In the following, (sin</> 9 , cos0 g ) will be abbreviated as 

( s g, c q)- 
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1. Self-energy of the top quark (see fig. [[]): 



<n}( 



H i 



■&*g* a C 2 (R) [2N{0,t,p)-4:m t F(0,t,p)} 

■^g*C 2 {R) [M(~g,t u p) ~ 2rh g s t c t F{g,t u p) 

+ M(g, t 2 ,p) + 2m g s t c t F(g, t 2 ,p) 
1 



16n 2 ' b t 



- sin 2 a 
2 

1 2 

+ -cos a 



N(H°,t,p)+m t F(H°,t,p)) 
N(h°,t,p)+m t F(h° J t jP ) 



+ - cos 2 [N(A, t, p) - m t F{A, t, p)) 

+ ^sin 2 /3[N(0,t,p) - m t F(0,t,p)} 

+ - [cos 2 /?M( J ff ± , b, p) + sin 2 /3M (0, 6, p) 



M(x°,ti,p)+m x o S<Ct (iV 2 +iV; 4 )F(x°,ti,p) 



| 2 M{x%t 2 ,p) - m x os t Ct(Nl + N( 4 )F(x Q v t 2 iP) 



h t UE [ c l I ^ T M(xt,h,p) + s 2 h I V i2 | 2 M(xtM,p) 



2. Self-energies of the stops (see fig. 0): 
(a) 

i<5 a ' 3 



in?, 



ii 



in 



ii 



aifi 



16vr 2 



16tt 2 



9tC 2 (R) 



2 

-p 2 - 2(m 2 + m?) 

8-KXfii, <, p) - 4G(#, t, p) - 8s t c t m t m g F(g, t, p)] 
1 



9tC 2 {R) 



qP 2 + - 4s 2 c 2 (m 2 - m 2 ) - 2K(0,t 1 ,p) 



+ L(0,*i,p) + (l-As 2 t c I t )K(t u t 1 ,0)+As 2 t c 2 t K(t 2 ,t 2 ,0) 
+h 2 \2(N + l> 2 c 2 (m 2 - m 2 2 ) + m 2 



*2 



16tt 2 
16vr 2 ' l t 



+ 2(JV + l^R^h, 0) + + ^ - 2Ns 2 c 2 )K(t 2 , h, 0)] } 



2/2-2 i 2 ~ 2 \ i 2 

^(c b m bi +s b m b2 ) + s t 



lK(b 1 ,b 1 ,Q) + s 2 b K(b 2 ,b 2 ,0)]} 

1 2 - 2 ,-'-2 2 ,-'-2 2/0,22 2/o 

-m H o sin a + -m h0 cos a + cos P + s t m H± cos P 



+ ^ sin 2 aK(#°, #°, 0) + - cos 2 aK{h°, h°, 0) 
+ ^ cos 2 /3i^(A, A, 0) + s 2 cos 2 PKiH*, H ± , 0) 
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16tt 2 



-2m 2 (< ) 2 F(#V~i,p) ~ 2m 2 (V 1 f) 2 F(H°,t 2 ,p) 
-2m 2 t {vSyF{h»M,p) - 2m 2 (V 1 h 2 °) 2 F(h°,t 2 ,p) 

-isin 2 /5(^) 2 F(At 2 ,p) - isin 2 /3(V;f ) 2 F(0,t 2)J 9) 



- sin 2 pp** ) 2 F(tf±, 6 1>P ) - sin 2 /?( V^f ) 2 F(tf ±, 6 2 , p) 

- sin 2 (3{vt ) 2 F(0, &!, p) - sin 2 flvg* ) 2 F(0, fc^ 

1 2 

9 — — ^ 



-p 
3^ 



£|iV j4 | 2 m 2 o 



-4]T I AT J4 | 2 K(x$,t,p) - 2£ I % I 2 G(xJ,*,p) 



16tt 2 



+ 2 StQ ]T(iV 2 4 + Nf 4 )m t m x oF(x° p t,p) 



h 



1 2 2 2 2 2 it/ 1 2 ~ 2 

-s t p - s t m b - s t }_^\ V i2 | m x + 



AslY, I ^ I 2 K(xt,b,p) - I ^ I 2 M 



30 9 2 S C 2 (R) 



i5 a P 
16tt 2 



16tt 2 
16tt 2 



-p 2 - 2(m 2 + m?) 



- 8if (<?, f , p) - 4G(g, t, p) + 8s t c t m t m g F(g, t, p)] 



<#7 2 (fl) 



1 2 , 1—2 i /i 2 2/ ~ 2 ~ 2 ^ 

gP + 2 m t2 + 4s t c ti m tl ~ m t 2 , 



2K(0,t 2 ,p) 



+ L(0,t 2 ,p) + (1 -4s 2 c 2 )ir(t 2 ,t 2 ,0) +4s^AT(ti,ti,0) 
+h 2 [-2(N + l)s 2 c 2 (m 2 tl - m 2 t2 ) + m\ 
+ 2(N + l)s 2 c 2 K(t 2 , t 2 , 0) + (4 + st- 2Ns 2 c 2 )K(t 1 ,t 1 , 0)] } 
h 2 {c 2 (c 2 b m 2 bl + s 2 b m 2 b2 ) + c 2 [c^C^, 6 1; 0) + s 2 b K(b 2 , b 2 , 0)] } 



1 2 - 2 . 1 2 2 . 1 2 2o,22 2 n 

-m H o sin a + -m h0 cos ct + cos P + c t m H± cos P 



+ ^ sin 2 aif(iJ , 0) + X - cos 2 aAT(fc°, h°, 0) 

+- cos 2 A, 0) + c 2 cos 2 PK{H±, if* 0) 

2 



-2mi{V£°) 2 F(H»M,p) - 2m 2 (V 2 ryF(H», t 2 ,p) 



H°\2 



-2ml{v£) 2 F{h»M,p) - 2m 2 t (V£) 2 F(h«,t 2 ,p) 
-isin 2 /3(^) 2 F(Ati,p) - isin 2 /?^ ) 2 F(0,t 1 ,p) 
- sin 2 (3(V 2 f ) 2 F(tf±, 6 1>P ) - sin 2 /3(F 2 f ± ) 2 F(tf ± , 6 2 , p) 
16 



h°\2 



ill: 



x° 

22 



22 



sin 2 /?(V 2 f ) 2 F(0,b 1 ,p)-sin 2 (3(V 2 T) 2 F(0,b 2 ,p) 

1 2 2 I at 1 2 ~ 2 



G±\2 ; 



16tt 2 "t 



-4]T I 7V j4 | 2 K(x°,t,p) - 2]T | 7V j4 | 2 G(*°,t,p) 
i j 

- 2s t c t J2(N] 4 + Nf 4 )m t m x oF{xlt,p) 



i5 a l 3 7,2 

16tT 2 



2 2 2 2 2 it/ 1 2 ~ 2 

-c t p -c t m b -c t }^\ V i2 | 



3 " 10 1 ^ 1 " 1 X; 
4c 2 £ I Vfc | 2 ^(x+,M-2c 2 £ | V a \ 2 G(xtAp) 



ill 9 



12 



in 



12 



in 6 12 



inf = 



l 12 



in* 



g 2 8 C 2 (R) cos20 t [-4m ( m 9 F(m 9 , 
|g {<#7 2 (i2) cos 20, - m 2 J 



16tt 2 



iS° 



16vr 2 



2stC t (K(t 1 ,t 1 ,0)-ir(t 2 ,t 2 ,0)) 



(N + l)s t c t cos20 t (m ti - m <2 ) 
+ (iV + l)s t c t cos 2(f> t (K(t 1: t 1: 0) - tf(f 2 , t 2 , 0))] } 
h t 2 {s t c t (c 2 m^ + s 2 m 2 2 ) 



16vr 2 "t 



+ s t c t 



h 



clKibubuO) + s 2 b K(b 2 ,b 2 ,0)]} 

m 2 H ± cos 2 (3s t c t + cos 2 /3s t c t K(H ± , if"* 1 , 0) 

^(^fVif )F(H°J llP ) - 2m 2 (V 1 fV 2 f)F(H°,t 2 ,p) 
-2m 2 {vEv£)F{h Q M,p) - 2m 2 (V^ 2 V 2 h 2 )F(h°,t 2 ,p) 

- sin 2 (3(V 1 H 1 ± V 1 H 2 ± )F(H ± , b u p) - sin 2 V 2 f )F(H±, b 2 ,p) 

- sin 2 p(V 1 fvg ± )F(0,b 1 ,p)-^n 2 p(V 1 fV 2 f)F(0,b 2 ,p) 



16vr 2 '"t 



h 



2Eft 2 - s 2 t N 2 A )m t m x oF(x° p t,p) 



ti 



'-s t c t p 2 - 2s t c t m\ - 2s t c t J2 I V i2 



2-2 

m + 

x7 



2s t c t X: I V l2 | 2 6,p) + G(x, + , M) 



In these expressions, we have used the following definitions for the Higgs-squark-squark 
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vertices in the (£1,^2) basis: 

V/( = sin cH s t c t (A t sin a — /xcosa) 

m t 

^12° = = o — cos2(j) t (A t sma — ficosa) 

VrS° = sin a s t c t (A t sin a — a cos a) 

m t 

V£ = Vi = {A t cotf3+u) 

Vg° = V 2 f = {A t -fi cot /3) 

= m t c t c b cot (3 + s t Cb(A t cot /3 + /i) 

V 12 = -m t c t s b cot /3 - s t s fe (v4 t cot (3 + y) 

= —m t s t c b cot (3 + c t c b (A t cot (3 + fi) 

V 22 = m t s t s b cot /3 - c t Sb(A t cot (3 + fi) 

V® = m t c t c b + s t c b (A t - (i cot /3) 

Vif* = -m t c t s b - s t s b (A t - y cot /?) 

= -m t s t c b + c t c b (A t - n cot (3) 

V 22 = m t s t s b - c t s b (A t - fi cot (3) 

The vertices can be obtained by ', replacing sin a — > cos a and cos a — > — sin a. 
3. Self-energy of the gluino (see fig. |]): 

• 2rafe 

m !W = -j^[2C 2 (G)N(0,~g,p)-4m g C 2 (G)F(0,~g,p)) 



^(p) = ~j^Y E [#(*i„*i,P) - 2m u s u c u F{t u ,U,p) + N(t 2i ,U,p) + 2m ti ^ Cti F(f 2j , * 
= "T^E [^(ftii.&i.P) -2m h s b ^F(b u ,b u p) + N(b 2i} b l} p) + 2m bi s h c bi F(b 2i 



In this expression, z is the generation index, and £, b stand for top-type and bottom- 
type. 

4. Self-energy of the Z° vector boson (see fig. (4|): 

n f z (p 2 ) 



n| (p 2 ) 

Km? z 



E ( ^"/^(/i. A. ^) + 2 (4 C /^ - s%c hR ) 2 m%F{h,k, Z) 
i=i lz 

+2(s 2 fi c fiL - c%c fiR fm%F{~f h , f\,Z) + ^\c\{c hL + c m ) 2 H(f n , f\, Z) 
~ 2 (44x + slc% R )mlF(f h J~ h ,0) - 2(44 £ + e^m^F^, f l2 ,0)} 





n|V) 

Km z 

nf(p 2 ) 

Km 2 z 



~H(H°, A,Z) + ± cos 2 (29 w )H(H ± , H±, Z) 



- \h{H\ H°, 0) - ±H{A, AO)- I cos 2 (29 w )H(H ± , F±, 0) 



In these expressions we have used: 

K 

CfR 



167T 2 



T 3 f — ef sin 2 9 
e f sin 2 9w 



w 



where / represents top- and bottom-type quarks, neutrinos and charged leptons (with 

This result 



e t = |,e 6 = -§,e„ = 0, e t = -1 and T 3t = T 3l/ = |,T 36 = T 3i = 
has already been obtained 0. We have reported it here for completeness, in the 
approximated form in which terms proportional to 0(m A z ) are neglected. 



Appendix B: Renormalization group evolution 

In this appendix we report the explicit solutions of the one-loop renormalization group 
equations for the parameters relevant to our calculations, in the approximation where the 
only non-vanishing dimensionless couplings are g s and h t . In the following, t = log(gj), where 
Qo is an arbitrary reference scale. The RGEs of interest can be extracted from the general 
formulae of [|14]. It must be stressed that, in accordance with our sign and phase convention 



0, the correct equations that must be solved correspond to the equations reported in the 
second paper of |14] with the substitution A q — > —A q . Following the notation of 0, we 



introduce a set of parameters: 



1 + ^a s (0)t\ (25) 



E = Z 3 9 (26) 

where b 3 = —3. It can be seen that, in our approximation, the integral F can be easily 
solved. Introducing the variable 

<t) = ^at(t) (28) 
the explicit solutions for a first set of equations are: 

a s (t) = a s (0)Z 3 (29) 
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a t {t) 


a t {0)E 




1 + ^(0)F 




fj. (t) 


= " (0) Uo)J 


1 

2 8 
r7 q 

^3 9 


v 2 2 (t) 


■ *»(S) 


1 

2 8 


Mt) 


= A t (0)(l-x)- 


-m 9 (0) 


rh g (t) 


= m g (0)Z 3 . 





16a a (0) „ l n , E 

— £Z 3 - -2 (i lk 

3 4vr 2 F ' 



The explicit solutions for the soft masses RGEs depend on the quantity: 

1 



Am 2 t = ^(m 2 Q (0)+m 2 T (0)+m 2 H2 (0))x 

1 „ / x , > r / e 

- -A t (0)x(l-x 



+ 



^m 2 g (0)x 



2 (t- - l) m g (0) - A«(0) 
£16a s (0) 



F 3 4tt 



tZ«-lt--l 



Then, we obtain: 





= m 2 Q {0)- 




Z 2 3 )m 2 g (0) 


m%{t) 


= m 2 T {0) - 




Z 2 3 )m 2 g (0) 


rh 2 B (t) 


= rh%(0)- 




Z!)m 2 g (0) 




= 1*4(0)- 


-Am 2 





3 s 

2 A 2 

-Am 2 

3 * 



(30) 
(31) 

(32) 

(33) 
(34) 



(35) 

(36) 

(37) 

(38) 
(39) 

We recall that the soft mass terms for the squarks of other generations (mQ i ,rhr i ,rhB i , with 
i — 1, 2) and (ffli^m^,^, with i — 1,2, 3) for the sleptons have been considered identical 
to mQ,m T ,mB, for simplicity. 
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